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We investigate multipartite entanglement for composite quantum systems in a pure state. Using the general- 
ized Bloch representation for n-qubit states, we express the condition that all fc-qubit reductions of the whole 
system are maximally mixed, reflecting maximum bipartite entanglement across all kvs.n — k bipartitions. As 
a special case, we examine the class of balanced pure states, which are constructed from a subset of the Pauli 
group V n that is isomorphic to Z£ . This makes a connection with the theory of quantum error-correcting codes 
and provides bounds on the largest allowed k for fixed n. In particular, the ratio k/n can be lower and up- 
per bounded in the asymptotic regime, implying that there must exist multipartite entangled states with at least 
fc = [0.189 7i J when n — > oo. We also analyze symmetric states as another natural class of states with high 
multipartite entanglement and prove that, surprisingly, they cannot have all maximally mixed fc-qubit reductions 
with fc > 1. Thus, measured through bipartite entanglement across all bipartitions, symmetric states cannot ex- 
hibit large entanglement. However, we show that the permutation symmetry only constrains some components 
of the generalized Bloch vector, so that very specific patterns in this vector may be allowed even though fc > 1 
is forbidden. This is illustrated numerically for a few symmetric states that maximize geometric entanglement, 
revealing some interesting structures. 



I. INTRODUCTION 

Quantum entanglement is certainly one of the most fasci- 
nating concepts arising in quantum mechanics, essentially be- 
cause it appears as a contradiction to reductionism, i.e., the 
principle by which understanding a complex system reduces 
to the description of each of its individual constituents fl~). 
As a matter of fact, a quantum composite system can pos- 
sibly be in a state such that its parts are more disordered - 
have a higher entropy - than the whole system. This peculiar 
property, known as the non-monotonicity of the von Neumann 
entropy J2j, is tightly linked to the notion of bipartite entan- 
glement. A pure bipartite entangled state, for example, admits 
a zero entropy, which translates the fact that one has complete 
knowledge about the joint system via its wavefunction. Its 
two parts, however, are mixed, so that they exhibit a nonzero 
entropy. In other words, one knows less about the parts than 
about the system taken as a whole, a property which cannot 
be conceived in classical terms. 

The essence of bipartite entanglement is thus that the in- 
formation about a quantum bipartite system is not only en- 
coded in its parts, but also in the correlations between them. 
Remarkably, when a bipartite quantum system is maximally 
entangled, the information appears to be fully encoded in 
these correlations and no longer in the system's constituents. 
Mathematically speaking, while the whole system is described 
as a pure state, its parts are individually described as maxi- 
mally mixed states (with a density matrix proportional to the 
identity). A paradigmatic example of such a situation is the 
Einstein-Poldolsky-Rosen (EPR) state of two qubits (3] |4], 
which is a pure bipartite state whose parts are maximally 
mixed: each qubit has an entropy of 1 bit, so its state is com- 
pletely unknown, while the 2-qubit joint state is perfectly de- 
termined. Equivalently, one observes that the entropy of one 
part conditionally on the other is negative (it is —1 bit), which 
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is a sufficient condition for bipartite entanglement and can be 
associated with a flow of (virtual) information backwards in 
time 0. 

A very intriguing question is of course whether similar sit- 
uations may exist if the system is made out of more than two 
parts. The underlying concept of multipartite entanglement 
becomes naturally much richer than bipartite entanglement, 
but also generally much more difficult to understand (see, e.g., 
[6 1 for a review on entanglement). It leads to stronger contra- 
dictions with local realism than bipartite entanglement]?), as 
well as to the existence of several inequivalent classes of en- 
tangled states even in the simplest case of three qubits |8). 
Multipartite entanglement is also crucial to applications, such 
as one-way quantum computing Q, and its dynamics when 
exposed to a dissipative environment has revealed a surpris- 
ingly large variety of flavors ifTOlfTTI . 

Among the possible approaches to multipartite entangle- 
ment, one of them consists of probing the presence of bipar- 
tite entanglement over all inequivalent bipartitions of all sizes 
ifTZll . Roughly speaking, the idea is to measure how much 
each subset of k out of n constituents (with < k < [n/2\) 
can be bipartite entangled with its n — k complementary con- 
stituents, knowing that there is a subtle balance with the bipar- 
tite entanglement exhibited by all other possible subsets with 
respect to their complements. This leads to the concept of a 
genuine multipartite entangled state, that is, a n-partite pure 
state such that none of its fc-partite subsets can be represented 
by a pure state (all subsets are mixed, hence bipartite entan- 
gled with their complements). 

One may even be more specific and seek for a strong form 
of a genuine multipartite entangled state. This would be a 
composite system in a pure state such that all of its individual 
constituents are maximally mixed, all pairs of its constituents 
are maximally mixed, all triplets of its constituents are maxi- 
mally mixed, and so on up to all fc-tuples of its constituents. 
This property, namely the fact of admitting maximally mixed 
reductions, is an ideal case of genuine multipartite entangle- 
ment. We expect that the constraint of having an overall pure 
state will set a limit on the possibility of having maximally 
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mixed fc -tuples for large k. Consider, for example, a system 
of n qubits in an overall pure state. It is tempting to search 
for n-qubit pure states having the property that all subsets of 
k qubits are maximally mixed up to size k = [n/2\. It has 
long been known, however, that except for a very few low- 
dimensional cases, such as the EPR state for n = 2, it is im- 
possible to find a n-qubit pure state exhibiting this strong form 
of genuine multipartite entanglement fPI . 

In this paper, we examine this fundamental question and 
investigate multipartite entanglement under two perspectives. 
First, we focus on the possible existence of n-qubit states that 
satisfy this property of admitting maximally-mixed fc-partite 
reductions for arbitrary n and k. We give examples of known 
states exhibiting this property, then provide a set of conditions 
that such states must satisfy as well as asymptotic existence 
bounds. The second part is centered on a weaker version of 
this maximally-mixed reduction property, which is motivated 
by recent results on symmetric states, a special class of states 
that are known to exhibit a high multipartite entanglement as 
measured in terms of their geometric entanglement lfT4l4T6ll . 

In Section II, we expose the representation of n-qubit pure 
states in terms of a generalized Bloch vector, which is very 
convenient in order to express the conditions that all fc-partite 
reductions are maximally mixed. This leads us to consider, 
in Sec. Ill, a class of n-partite pure states that we name 
"balanced". In the Bloch representation, they correspond to 
a Bloch vector with all components equal to a same value 
for indices belonging to some subset of the Pauli group (the 
other components being all taken equal to zero). This makes a 
connection with the theory of quantum error correction, from 
which we obtain lower and upper bounds on the highest al- 
lowed value of k for a given value of n. In particular, we 
show that there exist n-qubit states admitting all maximally- 
mixed fc-partite reductions with at least k = [0.189 nj when 
n — > oo. 

In Section IV, we then consider another natural class of 
states, namely symmetric states, among which it is known 
that some states with genuine multipartite entanglement can 
be found. We prove that, surprisingly, the symmetric states 
cannot have maximally-mixed reductions of size k that ex- 
ceed 1, regardless of n. In that sense, they are very far from 
the strong form of genuine multipartite entanglement that we 
seek. On the other hand, we show that the permutation sym- 
metry underlying symmetric states only puts constraints on 
the components of the Bloch vector with an even index, so 
that the components with an odd index may possibly be taken 
equal to zero. This brings us to investigate n-partite sym- 
metric pure states whose Bloch vector has many vanishing 
odd-index components. This investigation is carried out nu- 
merically, focusing on some symmetric states that maximize 
geometric entanglement as found in [15|. We show that some 
of these states are close to having maximally-mixed fc-partite 
reductions for large values of k (although, strictly speaking, 
k = 1), so that they approach the strong form of genuine mul- 
tipartite entanglement. Finally, some conclusions are drawn 
in Sec. V. 



II. MAXIMALLY-MIXED REDUCTION PROPERTY 

The property of admitting maximally-mixed reductions is 
encapsulated by the following definition: 

Definition 1. An n-qubit pure state \ip) is a k-MM state if all 
its reductions of size k are maximally mixed. Here and in what 
follows, MM stands for maximally mixed. 

Note that according to the definition, a fc-MM state is also a 
(fc — t)-MM state for < I < k, and, in particular, every pure 
state is a 0-MM state. A natural question which arises here is 
whether a fc-MM state exists for a given couple (n, fc). Note 
that no more than half of the qubits can be in a maximal mixed 
state, as a consequence of the Hilbert-Schmidt decomposition 
of the overall pure state. So, it is clear that fc-MM states cannot 
exist when fc > Here, we list the known facts about 

the existence of fc-MM states for small values of n: 



A. fc-MM states of small size n 

• For n = 2, it is easy to check that the four Bell states 

|4±) = M^n>, ,* ± > = ™>, ,i) 

are 1-MM states. 

• For n = 3, the Greenberger-Horne-Zeilinger (GHZ) 
state is a 1-MM state, while W state [ 8 ) is not. In gen- 
eral, for any size n, the generalized GHZ state 

\GHZ) . j^M , ,2) 

is a 1-MM state. 

• For n = 4, there exists no 2-MM state. To see this, 
consider the following 4-qubit states: 

\L) = ((1 - w)(|0011) + |1100)) + W 2 (|0101) + |0110) 

+|1001) + 1 1010) - |0000) - |1111))) , (3) 

\HS) = ^= (|0011) + |1100) + w(|0101) + 1 1010)) 
V6 

W(|0110) + |1001))) , uj = e 27ri/3 , (4) 

All of their 1-qubit reductions are maximally mixed, so 
they are 1-MM state, but this is not true for their 2-qubit 
reductions. \HS) was introduced in ifTTl and conjec- 
tured to be the 4-qubit maximally entangled state. Then, 
it was shown to be a local maximum of the averaged 2- 
qubit von Neumann entropy in |fl~8). Finally, in lfl9l . 
it was shown that the global maximum for all averaged 
2-qubit Tsallis and Renyi entropies is reached by \HS) 
for a < 2, and by the states \L) for a > 2. This implies 
that no 4-qubit state can have all of its 2-qubit reduc- 
tions maximally mixed, so indeed no 2-MM state of 4 
qubits exists. 
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• For n = 5, the two logical states |0)t and |l)x of the 5- 
qubit code introduced in [20 1 are both 2-MM state. It is 
easy to check that every qubits (or every pair of qubits) 
is found in a maximally-mixed state after tracing over 
the remaining qubits. 

• For n = 6, the four 6-qubit states constructed as logical 
Bell state using the previous 5-qubit code states, 



which, for a pure state, becomes simply 



1 



r a = —(i)\o a \^)- 



(9) 



If p is a quantum state, then hermiticity (p = p 1 ), positivity 



(p > 0), and normalization (Tr p 
constraints on the components r a 



= 1) give the three following 
ED: 



im^v _ |Q)|Q)l±|i)|i)l 

1 6 } V2 



\Mt 



|0)|l)L±|l)|0)i 
V2 



(5) 



r a eR,Va, 

r is in the positive cone 

r = i. 



are 3-MM states. 

• For n > 7, it is shown in [ 13 1 that no state can have all 
of its [n 1 2 J reduction maximally mixed. Note that the 
case n = 7 is not solved yet. 

We see that the 2-qubit 1-MM state and 6-qubit 3-MM state 
appear to be very special cases, with maximally mixed reduc- 
tions up to precisely half the number of qubits (the trivial up- 
per bound on k set by the Hilbert-Schmidt decomposition). 
In general, for an arbitrary number n of qubits, we can ex- 
pect that it will be possible to find a fc-MM state up to some 
threshold value /c max < L n /2J, which only depends on n. 



These three relations mean that, after translation by —1/2™ in 
the zeroth direction, a quantum state is completely represented 
by the vector r, which lives in the positive cone contained 
in a real subspace of a Schmidt-Hilbert space, the so-called 
generalized Bloch vector of dimension 4" — 1. Note that the 
concept of a positive cone embraces the idea that every convex 
combination of positive operators is also a positive operator. 

Pure states, i.e., rank-one projectors that satisfy p 2 = p, 
appear in this representation as vectors r that are constrained 
by 



a(3 



(10) 



B. Generalized Bloch vector formalism 

To analyze this question, we need first to introduce the gen- 
eralized Pauli matrices, which are the set of matrices con- 
structed in terms of all n-fold tensor products of the form 

0"a = °"ai ® °~a 2 ® ••• ® ^a„i (6) 

where each a a . represents respectively the 2x2 identity ma- 
trix or one of the usual Pauli matrices, depending on the index 
at E {0, 1, 2, 3}. The bold index a refers to a vector index 
containing the n indices on. There exist 4™ such matrices, all 
being traceless except for <r which corresponds to the 2™ x 2" 
identity matrix and admits a trace Tr(ero) = 2™. Using this 
set of generalized Pauli matrices {cr a }, it is also possible to 
construct the bigger set of the form {a a , —o- a ,icr a , —icr a }. 
This set of 4 n+1 elements becomes closed under matrix mul- 
tiplication and forms the so-called Pauli group V n - 

The set {a a } also forms a basis of a complex Schmidt- 
Hilbert space of dimension 4 n , so that every complex square 
2™ x 2™ matrix can be seen as a vector r in this space. For 
instance, a matrix p reads 



where g a p-f are the structure constants of SU(2 n ) defined as 
&a&f3 '■= 9ap~i °V Note that Eq. (jTOJ also automatically 
implies positivity. It can be decomposed in two independent 
relations 



M 2 = 

r) 



E^ 2 



2" - 1 
2 n + 1 



jk 



2™ 



2" 



(11) 
(12) 



where latin indices correspond to vector indices excluding the 
zeroth component {a = (0, i)), * is by definition the gener- 
alization of the cross product, and parentheses stand for sym- 
metrization, gujk) — (gijk + 5ifcj)/2. Relations (JTTJ and 
( fT2| ) express respectively that for a state being pure, its Bloch 
vector r should live on a sphere of radius R (which is actually 
the boundary of the positive cone) and should have a specific 
orientation. For more details on the generalized Bloch repre- 
sentation, see lETl . 



C. Conditions for maximally mixed reductions 



= r, 



while the components r a are given by the inverse formula 

1 



2" 



Tr(cr a p). 



(7) 

The generalized Bloch representation is a very useful tool 
in order to address the maximally-mixed reduction property. 
Let p be the density matrix of an n-qubit pure state \ip) living 
in the tensor product space H = C 2 <£> C 2 <E> ... <8> C 2 , and 
consider the bipartition H = Ha ® T~Lb where A and B are 

(8) 

defined as the sets of the first k qubits and last n — k qubits, 
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Pa (Computational basis) 



(Bloch representation) 



FIG. 1. The partial trace seen as a loss of information, both in the 
computational basis and in the generalized Bloch representation. 



respectively. In the generalized Bloch representation, the k- 
qubit reduced density matrix pa resulting from tracing out the 
qubits of B is given by 

p A = Tr B (p) = Tr B (^ r a <r a ) 

OL 

= ^2 ra A<XB a a A Tr B ((T aB ), (13) 

OLAO-B 

where the index a has been decomposed according to the bi- 
partion o.aolb- By using the fact that the matrices a aB are 
all traceless excepted for the one corresponding to the identity 
ohT-Ib, noted ao B , we get 

PA = ^2r aA o B a aA Tr B (cr 0B ) 
a a 

= 2 n - k J2r aA o B <? aA , (14) 
a A 

where r aA o B = r (aia2 ... Qfe00 ... 0) with n - k zeros in the 
vector index at the positions corresponding to B. If we now 
consider an arbitrary bipartition (A, B), we obtain an expres- 
sion similar to ( p"4| ) where in the vector index of r aA o B , the 
zeros are located at the positions of the traced out qubits. 
For instance, for a 6-qubit state, if the first, third, and last 
qubits are traced out, the component r aA o B corresponds to 

r (0a20a 4 a 6 0)- 

Equation ( fT4| i gives a very nice operational procedure for 
performing the partial trace. Usually, when expressing the 
density matrix in the computational basis, the entropy result- 
ing from tracing out a part of the system manifests itself both 
as a loss of some components and as the mixing of some other 
components of the original density matrix. In the generalized 
Bloch representation, the partial trace appears just as the loss 
of some components (see Fig.[T|. For instance, with the Bloch 
representation, the linear entropy of the reduced state simply 
reads 



S l {p A ) = l-Ti A {p A ) = l- 




Our concern now is the special case where all reduced den- 
sity matrices p A for a given size \A\ — k are proportional to 



the identity on Ha, '•<?., when the original state is a fc-MM 
state. In view of Eq. ( |14) , this is the case if all components r a 
with a vector index a containing at least n — k zeros vanish, 
except for the ro component which is always equal to 1/2™. 
By enumerating all these indices, we see that there are 




such components that must vanish. Let us define the weight 
Lu(a) of an index a (or by extension of a component r a or 
of a generalized Pauli matrix a a ) by its number of non-zero 
subindices. Then, we can establish the following proposition: 

Proposition 1. The n-qubit state p is a k-MM state if and 
only if its corresponding Bloch vectors r does not have any 
component r a with an index weight in the range < lo{ol) < 
k. 

This means that the vector r (after translation by —1/2™ in 
the zeroth direction) does not have any non-zero component 
in the subspace spanned by the basis vectors a a with index 
weight lower than or equal to k. If denote this subspace £). 
(of dimension Dk) and £k its orthogonal-complementary sub- 
space, we can establish the equivalent proposition: 

Proposition[TjDis. The n-qubit state p is a k-MM state if and 
only if its corresponding Bloch vectors r has zero components 
in £k, so that its support belongs to £/-■ 

The existence of a fc-MM state results from the compatibil- 
ity between having a pure state satisfying Eq. ( fT0) i and Prop. [T] 
at the same time. Such a compatibility is not easy to study for 
an arbitrary couple (n,k) without more information about the 
state. For this reason, we focus in the next Section on a class 
of balanced states that is suitable for analyzing this question 
and obtaining existence bounds. 

III. EXISTENCE BOUNDS FOR fc-MM STATES 

In the Bloch representation, constructing a pure state di- 
rectly in terms of the components of its Bloch vector involves 
the orientation relation ( fl2] >, which is difficult to manipulate 
in general. Even checking the purity of a given state numer- 
ically implies C(4 3n ) operations, and beyond 10 qubits the 
computation time becomes unreasonable on a standard desk- 
top computer. Instead, we will focus on a restricted class of 
states, which we call balanced states. 

A. Balanced pure states as fc-MM states 

Balanced state are defined in the Bloch representation as 
states whose non-zero components r a have all the same value, 
namely the same value as the identity component r . 

Definition 2. A n-qubit balanced state is expressed as 
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where S is a subset of the Pauli group V n , which entirely de- 
fines the state. 

Note that in this definition, the hermiticity of ps implies 
that S does not contain complex elements of the Pauli group 
of the form ±i a. We can now express the following theorem 
about the purity of such balanced states: 

Theorem 1. The n-qubit balanced state ps defined from the 
set S is pure if and only if S forms a group under matrix mul- 
tiplication that is isomorphic to ZJj. 

Proof. If ps is a pure state (p$ = p|) then we have 



2" ^ ' 



(18) 



This problem is fully equivalent to finding an additive self- 
orthogonal quantum error-correcting code over GF(4) 11231 . 
More generally, the connection between entanglement and 
quantum error correcting codes (QECC) was noted by several 
authors, and it was proven for example in lfl"3l that a QECC 
that can detect k errors is also a fc-MM state. The reciprocal 
of this statement can easily be understood by interpreting each 
Pauli matrix as an error operation. Indeed, for every tensor 
product of Pauli matrices with index weight < ixj(a) < k, 
Eq. (|9]l gives 



0. 



(21) 



for a fc-MM state \ijj), which means that each error a a is de- 
tectable because it rotates the state in an orthogonal subspace. 



a,r£S 



The uniqueness of the expansion ( fT7] i implies that ( fT8] > is sat- 
isfied only when S is closed under matrix multiplication, i.e., 
when S is a subgroup of V n . It follows that 



crGS TSSaeS 



(19) 



which implies that for ps being pure, the order of S should be 
2™. In summary, S must fulfill the three following properties: 

• S has a finite order equal to 2 n ; 

• S is abelian because normalization and hermiticity of 
Ps imply — er ^ S and ±i er ^ S, respectively; 

• All the elements of S have order 2, i.e., they are such 
that a 2 = do, since ±ia ^ S. 

According to the fundamental theorem of finite abelian groups 
El . this implies that S ~ Z 2 x Z 2 x ... x Z 2 = Z%. □ 

A simple example of such a pure balanced state is the state 
|0) = 1 00 • • • 0). Its components in the Bloch representation 
are given by Eq. dSl as 



r Q = i r Tr(a Q |0)(0|) = ^ r (0|a c 



— — (c a )oO) 



which are non-zero if and only if a a is a tensor product of 
identity and az matrices. There are 2" such matrices a a , 
which form the set S, and thus the state can be written explic- 
itly as 



|0> <0j = — (<7(00-00) +C(00---03) +O"(00-30) + 
' ' ' + f(33- -30) + (7(33. ..33)) j 



(20) 



where it is clear that S ~ Z2 just by relabelling the identity 
matrix <tq as "0" an d the Pauli matrix az = C3 as "1", for 
instance (T(oo— 03) — > (00 • • • 01). 

According to Prop. 1, a pure balanced state will be a fc- 
MM state if its group S does not contain elements with an in- 
dex weight lower than or equal to k (except for weight zero). 



B. Quantum Gilbert- Varshamov and quantum Hamming 
bounds on fc-MM states 

By exploiting this relationship, known bounds in the con- 
text of QECC can be mapped onto existence bounds for fc-MM 
states. The quantum Gilbert- Varshamov (GV) and quantum 
Hamming (H) bounds introduced in [24] give, respectively, 
lower and upper bounds on the number of errors that a quan- 
tum code can detect, for a given number of qubits. They can 
be directly translated in our context and give 



D [k/2 } + 1 < 2" < D k + 1. 

(H) (GV) 



(22) 



in terms of the subspace dimension D k . Asymptotically, these 
bounds can be written as [241 



(23) 



k \ <Q< f( k 
n I ~ ~ \2n 



(GV) 



(H) 



where 



f{x) = l-x log 2 3 + x log 2 a; + (l-x) log 2 (l-x), (24) 

is a decreasing function which has a root in xo ~ 0.18929. 
Physically, this means that it is always possible to find a n- 
qubit fc-MM state (n, k —> 00) that is such that, by keeping 
less than 19% or more than 81% of its qubits, we completely 
lose the information on the initial pure state. On the other 
hand, it is impossible to find such a state if we keep between 
38% and 62% of its qubits. The situation in the region be- 
tween 19% and 38% (or between 62% and 81%) is unknown, 
see Fig. 2. An intriguing physical implication of these bounds 
is that the entropy behaves as an extensive quantity (it is pro- 
portional to the number of qubits) in any subsystem as long 
as it has a size lower than 19% of the total system. It is only 
beyond this bound that, at some point, we observe a defect of 
extensitivity which originates from the purity of the state of 
the total system. 

Note that there exist more accurate bounds on QECC (see 
for instance (26)), but these are useless in our case. In- 
deed, finding a fc-MM state is equivalent to finding a one- 
codeword's code, and typically these other bounds become 
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stronger only when the number of codewords exceeds 1. In 
our case, there are no known better bounds than the quantum 
Hamming and quantum Gilbert- Varshamov bounds. 



C. Numerical bounds on fc-MM states 

Constructive upper and lower bounds can also be ob- 
tained numerically. For instance, the search of additive self- 
orthogonal quantum codes over GF(4) based on linear pro- 
graming has been performed up to around 100 qubits. Such re- 
sults are plotted in Fig. 2, based on Markus Grassl's database 




FIG. 2. (Color online) Constructive upper (thick red line) and lower 
(thick green line) bounds for the existence of fc-MM states based on 
Markus Grassl's database 1251 . The asymptotic limits on the domain 
of existence are also shown, namely the quantum Hamming bound 
(thin red line) and quantum Gilbert- Varshamov bound (thin green 
line). 



IV. SYMMETRIC fc-MM STATES AND BEYOND 

A. Symmetric states 

In 11271 . it is conjectured that any n-qubit state of maximal 
multipartite entanglement should be a 1-MM state, when the 
entanglement is measured through (the sum of) the negativity 
over all inequivalent bipartitions. We may naively extend this 
conjecture by saying that any n-qubit state of maximal multi- 
partite entanglement should be a fc-MM state, with k being the 
maximum allowed value as analyzed in the previous Section. 
While this is very well possible, we will see that the symmet- 
ric states, as defined below, are not the good candidates to test 
this conjecture. To understand why, we first recall that a sym- 
metric state is a n-qubit pure state that is invariant under any 
permutations of its qubits, that is 



U^) = |^),Vvre S n , 



(25) 



group of n objects. These symmetric states form a (n + 1)- 
dimensional symmetric subspace of H, which is often con- 
sidered as a good subspace to look for genuine multipartite 
entangled states, particularly in terms of their geometric en- 
tanglement lfl4l4T6l . Surprisingly, we observe that they satisfy 
the following theorem: 

Theorem 2. Symmetric n-qubit states are at most 1-MM 
states. 

Proof. We first express the unitary transformation that effects 
the transposition between qubits i and j, i.e., the SWAP op- 
erator. In terms of generalized Pauli matrices, this operator 
reads 



U, 




(26) 



In this expression we use the more usual notation for gener- 
alized Pauli matrices, i.e. the positions of non-identity matri- 
ces are indicated as superscripts while the positions of identity 
matrices are implicit, and the global identity is still written uq. 
For instance, a\\ corresponds to the action of a x on the first 
and second qubits. Now, we can write its expectation value in 
a given pure state as 




(27) 



Because the symmetric group S n can be generated by the set 
of all transpositions of two elements, the set of relations ( |2"5j ) 
is fully equivalent to the set of relations 

3 1 

L/yM = = ~> Vi > j € [l,n], (28) 

a=l 

which provide conditions on some components of the gener- 
alized Bloch vector r of a symmetric state At the same 
time, a fc-MM state with fc > 1 should have all components 
r ofc = with * 7^ h according to Prop. 1, which necessarily 
contradicts some of the relations (f28b. 



where is the unitary transformation that effects the per- 
mutation 7r over the set of qubits and S n is the symmetric 



□ 

Note that in the context of QECC, Theorem|2]means that it 
is not possible to find a code in the symmetric subspace detect- 
ing more than one error (or correcting any error). More phys- 
ically, we see that the permutation symmetry creates some 
frustration, which prevents the subsets of 2 qubits (or more) 
to be all maximally mixed. A natural question is of course 
whether permutation symmetry also manifests itself by con- 
straining higher index weight terms to be non-zero. To answer 
this question, let us define the concept of parity of an index a 
(and by extension of a matrix a a or of a component r a ): 

Definition 3. Let index a. contain Ai subindices 1, A2 
subindices 2, and A3 subindices 3. The parity of a is defined 
as even if the A; 's are all even. Otherwise, it is defined as odd. 
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For instance the indices (Oil) and (1122) are even, but 
(122), (0123) and (1123) are odd. Note that any component 
r a with an odd index weight oj (a) is necessarily an odd com- 
ponent. We can now state that 

Theorem 3. For a symmetric state, the even components r a 
of a given index weight w(a) cannot all vanish. 

Proof. Let us consider the expectation value of m transposi- 
tions acting on disjoined supports (m < [n/2j): 

m m / 3 \ 

fe=l k=l \ a k = l / 

-t / m 3 

\ fc=l os,=l 



(29) 



where, in this expression, there are (™) terms of the form 



E 



jii3i l '232---H3t 
aiaia2d2'"0't 



(30) 



which correspond to all the combinations of t transpositions 
chosen among the m transpositions. By induction, starting 
from the case t — 1 corresponding to Eq. ( |28] i, we conclude 
that all terms labeled by even indices should be equal to one, 
leading to the constraints 



E r hjii232—hjt — jL 
a 1 a 1 a 2 a2---ata t on.' 



(31) 



for t — 1, 2, ... m. Thus for symmetric states, even compo- 
nents r a of weight 2t cannot all vanish. □ 

Note that a priori the odd components are not constrained 
explicitly by the permutation symmetry. To see this, a reason- 
ing similar to the above proof can be done by considering the 
expectation values of a 3-cycle acting on qubits i\, i-i and 
We have 

3 3 

(U ili2 U i2i3 ) = ^\(a + J2<t)(<ro + ^4rm 

a=l 6=1 

(3 3 
i+E^kri^ + E^i^n^) 
a=l 6=1 

+ E^i<^r^>) 

ab ) 

/ 3 3 

= \ ( i + J2w<tw + X>kr iv) 



E^I CT - 3 ^)+ E £ <^<« 
a^6,e 



^2^3 I 
c6 



a—b 



where e Q b c stands for the completely antisymmetric symbol. 



Since we must have (Ui^Ui^) 
satisfied, it follows 



1 and since Eq. d28j) is 



E e " b 



E ^ ahcT ac 



1213 _ n 



(33) 



Since r* 1 ! 2 * 3 is completely symmetric in its lower indices for 



a symmetric state, Eq. (33 1 is necessarily satisfied. Then, we 



can generalize the last procedure by averaging the product of 
disjoined transpositions and a 3-cycle (acting on a disjoined 
support). Proceeding by induction in analogy to the reasoning 
leading to Eq. pT) , we obtain again a set of relations 



E 

a 1 a.2---a t 
616263 



e b 1 b 2 b 3 r : 



Jl3l«232---»tJtJlJ2»3 

aiaia2a2---a t a t bib2b-- 



= 



(34) 



(32) 



which are necessarily satisfied. Thus, in contrast with the even 
components which cannot all vanish in a symmetric state, the 
odd components are not constrained by permutation symme- 
try. 

It must be stressed that the odd components may, however, 
be possibly constrained for another reason, in particular as a 
result of the purity constraint ( fTO) . It is nevertheless natural to 
seek for symmetric states such that a large number of their odd 
components would vanish, which would correspond to high 
multipartite entanglement. In order to test this possibility, we 
now investigate specific states defined in lfl4T - fl6l . known to 
exhibit high geometric entanglement. 



B. Symmetric states with high geometric entanglement 

In EHED, some symmetric states with high geometric en- 
tanglement have been found. In particular, in |fl5l , an op- 
timization procedure was performed up to n = 12 in or- 
der to find the states that maximize their geometric entangle- 
ment. Those states are good candidate to test whether odd 
component r a indeed vanish. Of course, a direct calculation 
of these components by using expression |9| in the full 2 n - 
dimensional space is not realistic if n is not very small. For- 
tunately, by exploiting permutation symmetry, we obtain the 
following simplifications enabling an efficient calculation. 

The first simplification comes from counting the number 
of distinct components that we need to calculate for a sym- 
metric state. It is easy to see that U^pU\ = p implies 
r 7r(a) = r a for any permutation it £ S n . Thus, it is more 
convenient to label each component r a by a 4-component 
vector A = [A , Ai, A2, A3] which enumerates the numbers 
of subindices 0, 1, 2, and 3 in the index a. For instance, 
we have r( ii 3 ) = r (1031) = r[ 1>2)0 ,i]. For a given A, the 
number of equal components r a corresponding to the same 

is given by the multinomial coefficient f. . . ) = (j) . 
Thus, instead of having to calculate 4™ components, only 
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(CD) = rr 1 ) = > + 3)(» + 2)(» + 1) * £ dis- 
tinct components are needed, corresponding to the number of 
distinct multinomial coefficients. Note that permutation sym- 
metry also implies that U w pUL = p for any permutations 
7T, 7r' € £„, which leads to other constraints on the compo- 
nents r a . 

The second simplification comes from calculating Q in the 
so-called Dicke basis, which is a natural basis of the (n + 1)- 
dimensional symmetric subspace and allow us to use the states 
as expressed in this basis in ifTBll . Any symmetric state can be 
decomposed as 



ft=0 

where \S 7 k l ) are the Dicke states, 4 6C and Ylk=a \dk\ 2 = 1. 
The Dicke states are written in the computational basis as 



\st) 



-1/2 



E 

7TSS„ 



C/J00- 



-1/2 



(35) 



where j is a binary vector of size n such that | j | = X)"=i 3i — 
k. Symmetric states can also be defined in terms of the Ma- 
jorana representation [28]. In this representation, every sym- 
metric state is characterized by a collection of n one-qubit 
states \qi) — ccj|0) + which can be viewed as n points in 
the surface of the Bloch sphere, according to the expression 



e i0 

7TGS„ 



Un\qi)\q2) ■■■Wn 



(36) 



for some phase 9 and some normalization factor Af. We can 
move from the Majorana representation to the Dicke basis 
thanks to the relation |29l 

/ \ -1/2 

^ k = ( h ) E ' ' ' y^(k) x ^(k+l) ■ ■ ■ »7r(n)) 

but this is really inefficient as it involves a sum over all per- 
mutations. Instead, we use the fact, also noted in |[T4l [29), 
that the Majorana parameters Zi = Xijyi are the roots of the 
polynomial 



F(z) = £(-!)< 



k=0 



-1/2 



(37) 



Thus, knowing the Majorana parameters Zi, the Dicke com- 
ponents dk can be calculated in about 0(n 3 ) operations by 
solving the set of linear equations J2k Aikdk — 0, with A 

being a matrix of entries An, = (— l) fc (?) 2*. 

The last simplification concerns the expression of the gen- 
eralized Pauli matrices themselves. Indeed, in order to com- 
pute the components rv of a symmetric state of known com- 
ponents dk in the Dicke basis, we only need the symmetric 



part of the generalized Pauli matrices a a , that is, their projec- 
tion into the symmetric subspace. These symmetric matrices 
are (ra + 1) x (ra + 1) matrices in the Dicke basis, which we 
note as 7v. Just as for the components r a of a symmetric 
state, there are (?) generalized Pauli matrices a a which are 
all projected onto the same symmetric matrix rr, labeled by 

the index A. The matrix elements of rr in the Dicke basis are 



1/2 / x -1/2 



E E <j>«ij"> 

lil=& |j'i=fc' 



y> (j>® A ° of Al g a? A2 af X3 \j') 

lil=* V (fe) 

li'l=* 



where we have taken an arbitrary order for the individual Pauli 
matrices in cr a . Then, the vectors j and j can be cut in four 
pieces j a and j a of size A Q (a = [0, 3]). To simplify the nota- 
tion, we note as J the domain satisfying the set of constraints 
on vectors j a and j' a . It follows 



( t a w = E 



(joK® A °Uo)Oik f Al \3i) ■ ■ ■ <J 3 k 3 8A 'l Ja) 
.fe) (fc') 



.7 



By using the definition of the usual 2x2 Pauli matrices, we 
can write each factor as 



(iikf Al b"i) = ^ 1 jr J 

(32Wf X2 \j'2) = (-l)^W 2 S j2 jr, 
(3 3 Wi X3 \33) = (-l) lj > l S j3 f a , 



(38) 



where barred vectors stand for the complementary vectors, for 
instance (1011) = (0100). These equations can be rewritten 
in terms of the new indices k a = \j a \ and k' a = \j' a \ as 



{k l \vf Xl \k' 1 )=6 kl (Al - fe;) , 

(* 2 k? Aa l^) = (-i)^W Aa ^(A a -*i), 

{ (h\*i x *\k' 3 ) = (-i) k, *5 k3k ,, 



(39) 



and the sum over each j a can be replaced by a sum over each 
k a weighted by a factor ( Aq ) . Eventually, the matrix elements 
can be reexpressed as 



fa) 



A A' 



E 



) -(2fc 3 +3A 2 -2A 2 ) Ao\ /AA fA 2 \ (\ 3 



Mo\ /AA /A 2 \ /A 3 \ 
IaV Vfei/ Vfc 2 / Vfc 3 / 



\ n, / \ n, / 



(40) 



where the sum is taken over the four indices ko, ki, k-i and 
kj, which can take values between and Ao, Ai, A2 and A3, 
respectively, with the two constraints 



k = k + k\ + k 2 + fc 3 , 
k' = k + (Ai - ki) + (A 2 



(41) 
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Note that in the worst case where each X a ~ / 4J , this cal- 
culation implies to calculate about 0{n A ) terms. As an exam- 
ple for two qubits (n = 2), the symmetric part of the general- 
ized Pauli matrices 0713) and crpx) corresponds to 



/ o 



r [0, 1,0,1] 



l 




1 


1 

V2 



(42) 



By using explicit expressions of symmetric states with max- 
imum geometric entanglement taken from ifTBI or 3D coordi- 
nates available in the Sloane database IT301 . these simplifica- 
tions allow us to calculate efficiently the distinct components 
in the Dicke basis as 



1 



For n = 4-12 (and also for n — 20), the ratios between the 
number of zero odd components and the total number of odd 
components are gathered in Table [I] We observe that a large 



n 


State 


Zero odd/Total odd 


4 
5 
6 
7 
8 
9 
10 
12 
20 


|^4) = ^|5o 4 ) + V / ||5 3 4 ) 

|Vs) ~ 0.547|S^) + 0.837|S|> 
\^ 6 ) = (\St) + \Sl))^2 
\4, 7 ) = (\S 7 1 ) + \Sl))V2 
|Vs) ~ 0.672|S! 8 ) +0.741|Sf) 
|^ 9 ) = (|5 2 9 ) + |S?))V2 
IV-io) = (|S 2 10 ) + |S 8 10 »V2 
IV12) = (|5 2 10 ) + |S 8 10 »V2 
Dodecahedron state from BOI 


18/25=72% 
36/46~ 78 % 
64/64=100 % 
90/100=90 % 
94/1 30~ 72% 
164/185~ 89% 
230/230=100% 
34 1/37 1~ 94% 
1266/1484 ~ 85% 



TABLE I. Proportion of vanishing odd components in symmetric 
states that maximize geometric entanglement (from 1 15 ]). Note that 
\tpio) is not the maximum but only a state really close to it which 
allows an explicit writing. 

proportion of odd components vanish for these states. In the 
special cases n = 6 and n = 10, really all odd components 
vanish. For the case n = 4 and n = 12, the states con- 
structed thanks to the 3D coordinates available in [ 30 1 give 
better results, namely the ratio 24/25 ~ 96% for n = 4 and 
371/371 = 100% for n = 12, even though these states are 
equivalent to those of [ 15 1 (i.e. related by symmetric unitary 
transformation £/®"). This is the case because the proportion 
of vanishing odd components is basis-dependent, in the sense 
that two equivalent symmetric states will have in general a 
different structure in the vector even if they have the same 
entanglement content. 

Finally, we observe that the structure of the vector r? often 
takes a particularly simple form, especially for these states 
with a large proportion of vanishing odd components. For 
n = 4 and n = 6, all the non-zero components are given by 



r 7r[4, 0,0,0] 
r 7r[2,2,0.0] 
' r n[l, 1,1,1] 

' r 7r[6, 0,0,0] ' 
^[4,2,0,0] 
. r 7r[2,2,2,0], 




(43) 



(44) 



In these expressions, the permutation symbols ir applied on 
vectors A indicate that components with the same vector in- 
dex A up to some permutation are equal (or of opposite sign). 
For instance, for the 4-qubit state, it means that rvoon) = 
r (oo22) = ~ r '(22ii)- Note that the presence of the minus sign 
can be see as a consequence of the purity constraint (jT0j». For 
n = 10 and n = 12, the structure is similar even if a little bit 
more complex because it involves one or two different values 
per permutation of the A index. In particular, for n = 12, the 
structure is similar to that of n = 4 or 6 in the sense that we 
have the non-zero components: 



n = 12 



/ r Tr[12,0,0,0] N 
r rr[10,2,0,0] 
r 7r[8,4,0,0] 
V r 7r[8,2,2,0] 



/ 1 

±1/3 
2/10 
V±l/15y 



(45) 



We believe that a possible new approach to analyzing max- 
imum entangled states in the symmetric subspace for higher 
values of n should be inspired by these nice structures, and 
take the vector as a starting point. 



V. CONCLUSION 

We have investigated the maximally mixed reduction prop- 
erty through the concept of fc-MM states. By making use 
of the generalized Bloch representation in which the partial 
trace operation takes a simple form, we expressed the condi- 
tion that a fc-MM state must satisfy in terms of its generalized 
Bloch vector components. Considering the class of balanced 
fc-MM states and a connection with quantum error-correcting 
codes, we found asymptotic lower and upper bounds on the 
reduction size k for n — > 00. Then, we analyzed the class 
of symmetric states, which led us to consider a weaker ver- 
sion of the maximally mixed reduction property. We showed 
that symmetric states cannot have maximally-mixed fc-qubit 
reductions with fc > 1, which is linked to the fact that some 
weight-two component of their generalized Bloch vector must 
necessarily be non-zero. In other words, symmetric states do 
not obey the maximally mixed reduction property (they cannot 
be fc-MM states with fc growing linearly in n). However, we 
showed that the constraint of admitting non-zero components 
only holds for even components, so odd components are not 
constrained by permutation symmetry. We studied the case 
of symmetric states which maximize geometric entanglement 
(up to n = 20) as examples of states admitting many zero 
odd components, witnessing a high multipartite entanglement 
content even though they do not obey the maximally mixed 
reduction property for fc > 1, This led us to observe some 
interesting structures in the Bloch vector of states maximizing 
the geometric entanglement, which may open new perspec- 
tives in the analysis of multipartite entanglement. 
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